The Origin of Neutral Hydrogen Clouds in Nearby Galaxy 
Groups: Exploring the Range Of Galaxy Interactions 
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ABSTRACT 

We combine high resolution N-body simulations with deep observations of neutral hydrogen 
(HI) in nearby galaxy groups in order to explore two well-known theories of HI cloud formation: 
HI stripping by galaxy interactions and dark matter minihalos with embedded HI gas. This 
paper presents new data from three galaxy groups, Canes Venatici I, NGC 672, and NGC 45, 
and assembles data from our previous galaxy group campaign to generate a rich HI cloud archive 
to compare to our simulated data. 

We find no HI clouds in the Canes Venatici I, NGC 672, or NGC 45 galaxy groups. We 
conclude that HI clouds in our detection space are most likely to be generated through recent, 
strong galaxy interactions. We find no evidence of HI clouds associated with dark matter halos 
above Mhi ~ 10 6 M Q , within ± 700 km s _1 of galaxies, and within 50 kpc projected distance 
of galaxies. 



Subject headings: galaxies: evolution, 
galaxies, cosmology: dark matter 

Introduction 



galaxies: interactions, ISM: clouds, ISM: evolution, radio lines: 



Neutral hydrogen observations of the Milky 
Way reveal a repository of cold gas clouds that 
do not appear to rotate with the galaxy. Since 
they are kinematically decoupled from galactic ro- 
tation, they are known as High Velocity Clouds 
(HVCs). Generally, a cloud is categorized as a 
HVC if its velocity deviates by more than 50 km 
s _1 from the range allowed by a simple model of 



differential galactic rotation (Wakker & van Wo- 
erden||1991[ ). 

Milky Way HVCs were first discovered by 



Muller et al. (1963), after detection of interstellar 



clouds through absorption line studies prompted 
systematic searches for HI clouds in equilibrium 



with a Galactic corona (Wakker & van Woerden 



1997). Subsequent studies have revealed that the 



demographics of HVCs in the Milky Way are rich 
and varied: HVCs are both ionized and neutral, 
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occur in long streams and spheroidal clouds, and 
are found in large complexes and as small isolated 
clouds. The neutral hydrogen HVCs cover nearly 



40% of the sky ( 


Wakker & van Woerden 


1997 


Lockman et al. 


2002 


). In contrast, the ionized 



clouds cover as much as 85% of the sky (Shull 



et al. 20091. In total, these clouds make up as 
much as 10% of the total HI mass of the Milky 
Way fWakker fc van Woerden|l997 |. 



Possible analogs to the Milky Way clouds have 
been found around nearby galaxies including M31 



(Thilker et al. 20041, M33 (|Grossi et al. 2008), and 


others ( 


van der Hulst fc Sancisi||1988 


Miller et al. 


2009 


Hess et al.|2009 


). The best-studied external 



HI cloud system in an isolated galaxy is that of 
M31. Around the galaxy, [Westmeier et al.| (|2008| 



found a population of HI clouds with properties 
consistent with those of Milky Way HVCs. The 
M31 clouds have typical HI masses of 10 5 M Q , and 
the system extends 50 kpc from the M31 disk. 

There are several proposed theories for the 
source of HI clouds. Broadly, they consist of two 
classes: those that invoke a galactic plane origin 
and rely on secular dynamics to remove the gas 
from the disk, and those that invoke an extra- 
planar origin and tie the dynamics to large scale 
structure. Table [T] lists the expected properties 
of HI clouds in each case, which we describe be- 
low. Note that the properties of cold-accretion 
generated clouds are not well-constrained by cur- 
rent simulations, particularly in terms of their ve- 
locities and sizes. 

In a galactic fountain, gas is blown out of the 
disk by a supernova outburst, after which it cools, 



condenses and rains back onto the disk ( Shapiro & 
Field|1976| pre gman|1980| . The galactic fountain 
model for HVCs is supported by the existence of 
a hot (10 6 K) galactic corona with a scale height 
of a few kpc, detected in x-rays around galaxies 
(e.g., Strickland et al. (2004)). This corona is sup- 



plied with gas by supcrnovae explosions from the 
disk, and radiative cooling causes the formation of 



HI clouds through thermal instabilities (Shapiro 



& Field 1976 Bregman 19801. Although galac- 



tic fountains most certainly exist at some level 
in all star-forming galaxies, not all observed HI 



clouds are produced this way ( Fraternali & Binney 



2006). In particular, the galactic fountain model 
fails to explain HI clouds further than a few kpc 
from the disk of a galaxy, beyond the extent of 



the hot galactic corona. Our observations are not 
sensitive to galactic fountain-generated HI clouds, 
as we are limited by both angular resolution and 
mass sensitivity. 

Galaxy interactions may generate HI clouds 
through tidal stripping (Putman et al. 2003 



Chynoweth et al.|2008| |Sancisi et al.|2008[ ). In the 



Milky Way, the Magellanic Stream is a prime ex- 
ample of this phenomenon-here gas was stripped 
out of a satellite galaxy during a minor merger. 
Another possible source of HI clouds far from 
galaxy disks is leftover fragments of tidal tails. 

HI clouds may be tracers of a population of 
10 7 — 1O 1O M dark matter halos (Kravtsov et al 



2004). Cosmological simulations based upon the 
accepted paradigm of ACDM (Cold Dark Matter) 
cosmology, such as the Via Lactea II simulation of 



Diemand et al. ( 2008 1 , consistently predict a high 
degree of dark matter substructure. This aspect 
of simulations is not well-matched to observations- 
for instance, the Milky Way is observed to have a 
factor of at least 4 fewer satellites (in the form of 
dwarf galaxies) than predicted by simulations, de- 
spite sophisticated treatment of baryonic heating, 



cooling, and feedback in low mass halos (Simon & 



Geha 2007). One possible solution to this "miss- 
ing satellites problem" could be HI clouds, if such 
objects are tracers of dark matter halos. How- 
ever, there are still not enough clouds observed to 
supply the predicted number of satellites, and the 
clouds that are observed do not have the expected 
spatial or kinematic distribution of dark matter 



halos (Chynoweth et al. 2009 2011). 



Numerical simulations show that gas may ar- 
rive to galaxies via accretion along dark matter 
filaments ( Keres et al.|2005 |. Gas accretion along 
filaments appears to happen in two modes, a "hot" 
(T ~ 10 6 K) and a "cold" (T < 10 5 K) mode. The 
cold mode of accretion, where gas is not shock 
heated to the halo virial temperature as it falls in, 
is of primary relevance to HI studies, since some 



cold-mode accretion gas may appear as HI (Keres 
fc Hernquist|2009 ). HI clouds could be a result of 



instabilities in a hot halo that causes the gas to 
cool in clumps and fall toward the galactic mid- 



plane ( Mailer & Bullock|2004 Keres & Hernquist 



2009). In this case, HI clouds are cold, pressure- 
confined clouds without an associated dark mat- 
ter halo (Keres & Hernquist 2009). When cold 
gas encounters the ambient hot halo gas in a 
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galaxy, a density inversion develops which causes 
a Rayleigh- Taylor instability. This scenario is sup- 
ported by the kinematics of the Milky Way HVC 
system, which imply significant inflow (Mirabel 



fc Morras|[l984l |Simon et al.||2006[ |Wakker et al. 
20071) 



Determining the relative importance of each HI 
cloud formation channel requires us to know the 
phase space and mass distribution of the clouds, 
and measurement of the physical properties of HI 
clouds is dependent on an accurate value for dis- 
tance. For example, mass scales as D 2 , size as D 1 , 
and density and pressure as D , where D rep- 



resents distance (Wakker & van Woerden 19971 



Unfortunately, the distances to the clouds in the 
Milky Way are very difficult to measure. HVCs 
are starless, so none of the traditional methods for 
measuring the distance to a star can be used. The 
clouds' anomalous velocities and extraplanar po- 
sitions make it impossible to use a model of disk 
rotation to determine their distances. The best 
that can be done in the Milky Way is to use stel- 
lar absorption lines to determine upper and lower 



distance limits for a cloud ( 


Schwarz et al. 


1995 


Wakker & van Woerden||1997 


). In order to deter- 



mine the upper distance limit, absorption in the 
spectrum of a star at the velocity of the cloud 
must be detected in a star in the line-of-sight of 
the cloud. For a lower distance limit, a signifi- 
cant non-detection of the absorption line in a star 
in the line-of-sight is required. Many attempts 
have been made to measure the distances of clouds 
in the Milky Way using stellar absorption (see, 



e-g-, 


Schwarz et al. 


( 


1995 




van Woerden et al. 


(1999 


); 


Ben Bekhti et al. 


2006)). These studies 



have placed cloud distances in the range of 5-10 
kpc above and below the disk ( |Barentine et al. 
2008). Because background stars are required for 



this method of distance measurement, only clouds 
that are relatively close to the disk can be mea- 
sured. 

An advantage to studying HI clouds in exter- 
nal galaxies is that there are many robust meth- 
ods for measuring the distance to a galaxy, so 
the distance to any associated HI cloud is known 
much more accurately. Additionally, the projected 
distance between a cloud and galaxy is easy to 
measure for extragalactic HI clouds. Galaxies are 
most commonly found in groups (e.g. Geller & 
Huchra| fll"983l; fFu\ty\ (fl987l); |Eke et al.| ([20041 



Tago et al. (2008)), so groups are logical envi- 
ronments in which to search for HI clouds. One 
HI cloud search of galaxy groups has been done 



- Pisano et al. ( 2007 ) searched a sample of loose 
groups that are Local Group analogs. They did 
not find any HI clouds, but used their observa- 
tions to set the constraint that HI clouds should 
be on average clustered within 90 kpc of group 
galaxies, with an average mass of 4 x 10 5 M . 



In Chynoweth et al. (20081, we presented ob- 



servations of the M81 group. In |Chynoweth et~aT 
(120091), we investigated the NGC 2403 group. In 



Chynoweth et al. (2011), we expanded the obser- 



vations of Chynoweth et al. ( 2008 1 and Chynoweth 



et al. (2009) to include the the filament between 
M81 and NGC 2403. In this paper, we search three 
more groups and synthesize this data in the con- 
text of our previous results. 

This paper is organized as follows. In Sec- 
tion [2j we describe our sample of galaxy groups 
and method of selection, introduce our method of 
quantifying galaxy interaction strength. In Sec- 
tion [3j we describe the observation strategy and 
data reduction process. Section [4] presents a dis- 
cussion on the completeness of the study. Section 
[5] describes the predictions of numerical simula- 
tions for HI clouds associated with dark matter 
halos. Finally, in Section [6] we present our results 
and conclusions, and discuss future work. 

2. Sample Selection 

We observed five nearby galaxy groups with 
the 100m Robert C. Byrd Green Bank Telescope 
(GBT) at the NRACQ in Green Bank, West Vir- 
ginia. We observed a wide angular area and veloc- 
ity range around each galaxy, and the properties of 
detected HI clouds were catalogued and compared 
against simulations. 

The five galaxy groups were chosen with the 
following criteria. 

1. Galaxy groups must span the range in inten- 
sity of galaxy interactions. 

2. Galaxy groups must be nearby, for adequate 
angular resolution and sensitivity. We aimed 



lr The National Radio Astronomy Observatory (NRAO) is a 
facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 
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for a maximum group distance of 7.5 Mpc, 
corresponding to an angular resolution of 20 
kpc. 

3. Galaxy groups must be at least partially dis- 
tinct from the Milky Way in velocity space 
to avoid confusion with foreground HI gas. 

Criterion 1 is essential in order to explore the 
role of interactions in the generation of HI clouds. 
We attempted to choose galaxy groups with dif- 
ferent levels of galaxy interactions. Although in 
some cases it is obvious that a galaxy interaction 
is taking place, it is not always a simple mat- 
ter to determine whether a particular galaxy has 



undergone a recent interaction (Conselice et al 



2008 Darg et al. 2010a). Determining the rela- 



tive strength or intensity of interactions is even 
more difficult. However, we have made a first 
attempt towards quantifying galaxy interaction 
strength within galaxy groups by creating "in- 
teraction indices", which we have applied to our 
galaxy groups in order to rank them in order of 



interaction strength (see Section 2.6 1. 

Criteria 2 and 3 are difficult to reconcile, so the 
two closest groups (M81 and NGC 2403) partially 
overlap with the Milky Way in velocity space. The 
groups discussed in this paper do not overlap in 
velocity with the Milky Way. The galaxy groups 
observed for this study are described below. Table 
2] lists all galaxy groups within 7.5 Mpc (iFouquel 



et al. 1992 Tully 1987). Our observations cover 



63% of galaxy groups within 7.5 Mpc. Figure [T] 
shows a diagram of the distribution of the observed 
groups. Here we list the groups. 

2.1. M81 Group 

The M81 group is one of the closest galaxy 
groups to the Milky Way at a distance of 3.5 Mpc. 
It is comprised of 4 major galaxies (M81, M82, 
NGC 3077, NGC 2976) and over 40 dwarf galaxies. 
Although the galaxies are optically undisturbed, 
HI observations reveal that the galaxies are un- 
dergoing strong interactions. Additionally, M82 
is an extreme starburst galaxy and many of the 
dwarf galaxies show morphological and kinematic 
irregularities. Table [3] lists galaxy properties. 



2.2. NGC 2403 Group 

The NGC 2403 group is located at approxi- 
mately the same distance as the M81 group. The 
NGC 2403 group contains the bright spiral galaxy 
NGC 2403, 3 additional large galaxies (NGC 2366, 
UGC 4305, UGC 4483), and at least 4 dwarf or 
low surface brightness satellite galaxies. The NGC 
2403 group does not show clear signs of interac- 
tion. However, some of the group galaxies show 
subtle irregularities which may be due to previous 
interactions between group galaxies. The largest 
galaxy, NGC 2403, shows a thick, lagging HI layer 
and a few HI filaments with anomalous velocities 
that are similar to the M81/M82 clouds in mass 
(Fraternali et al. 2002[ ). These features may be 



analogues to Milky Way HVCs, and their origins 
are similarly unclear. They have been attributed 
to a combination of outflow from a galactic foun- 
tain and accretion from the intergalactic medium 



(Fraternali & Binney 2006); both of these phe- 



nomena may be induced by interactions between 



group galaxies (Larson & Tinsley 1978 Sancisi 



et al. 


2008 


. UGC 4305 has a come 


tary appear- 
), which may 


ance in HI 


Bureau & Carignan||2002 



be caused by ram pressure from the intragroup 
medium. The dwarf galaxy DDO 53 is also kine- 
matically irregular in HI ( Begum et al.|2006 1 . The 
other galaxies in the group appear to be normal. 
Table [4] lists galaxy properties. 

2.3. Canes Venatici I Group 

The Canes Venatici I (CVn I) galaxy group is 
located approximately 4 Mpc away and contains 
41 galaxies distributed in a diffuse cloud. The 
brightest galaxy, M94, is apparently isolated but 



is morphologically disturbed and has a ring ( Mul 
der fc van Driel||l993 ). Five of the galaxies show 
bursts of star formation (Kaisin & Karachentsev 
2008). Table [5] lists galaxy properties. 



2.4. NGC 672 Group 

The NGC 672 and NGC 784 groups form a 
6° linear filament of galaxies at a distance of ap- 



proximately 5 Mpc ( |Zitrin fc Brosch||2008[ ). The 
galaxies in both groups show evidence that a burst 
of star formation occurred approximately 10 Myr 
ago. Zitrin & Brosch (2008) conclude that these 



coincident SF bursts are most likely not due to 
galaxy interactions, since the SFR of the galaxies 
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does not correlate with projected distance from 
the brightest and most massive galaxy. Instead, 
they attribute the bursts in star formation to ac- 
cretion along the filament. However, it is well 
known that NGC 672 a nd IC 1727 are undergo- 
ing a strong interaction ( Combes et al.|l 9 80 ; Sohn 
& Davidgc 1996). Table 6] lists galaxy properties. 



2.5. NGC 45 Group 

The NGC 45 group contains 3 large, low sur- 
face brightness galaxies (NGC 45, NGC 24, NGC 
59) and is located at a distance of approximately 
6 Mpc dFouque et al.||1992|). NGC 45 has a kine- 



matical twist in the major axis, but NGC 24 is 
kinematically regular ( Chemin et al.|[2006 |. How- 
ever, NGC 24 has star forming regions scattered 
across the entire disk (Rossa & Dettmar 2003| ). 
NGC 59 is the third galaxy in the group, and it 
appears to be normal and undisturbed. Table [7] 
lists galaxy properties. 

2.6. Quantifying Galaxy Interactions 

Measuring the probability of a recent galaxy 
interaction or merger in an observational dataset 
has been a goal for decades; common methods em- 
ploy pair counts ( Zepf fc Koo|1989 Carlberg et al 
T9Ml|Patton et al .|1997||Le Fevre et al.|2000||P~ 



ton et al.||2002l |Conselice et al.||2003"l |Bell et al 



2006 Berrier et al. 20061, morphological classifi 



cation techniques (i.e., concentration, asymmetry, 



and dumpiness (Abraham et al. 1996 Conselice 
et al.| 20031 or the Gini coefficient (Lotz et al. 



2004)), or star formation (Sanders et al. 1988 



Knapen fc James||2009[ ). None of these are iron- 
clad merger indicators- phenomena that are some- 
times induced by interactions are not a guarantee 
of interactions. Furthermore, it is not clear how 
much weight to give each possible indicator. With 
these caveats, we used a combination of these pa- 
rameters as a guide to rank the level of interac- 
tions in a galaxy group. Here we describe how 
we compiled merger characteristics into an "inter- 
action index". We use this as a rough guide to 
quantifying galaxy interactions. 

A key symptom of galaxy interaction is found 
in the galaxy morphology. Various morphological 
features may be induced by galaxy interactions. 
Although spiral, barred, and elliptical galaxies can 
all be related to interactions, the most reliable 



morphological indication of recent interactions are 
tidal bridges, tails, and shells. 

Secondly, strong central star formation is a key 
feature of recent or ongoing interactions. Many in- 
teracting systems show indications of starburst ac- 
tivity such as high levels of Ha emission and high 
infrared luminosity (e.g., Sanders et al. (1988); 



Knapen & James (20091). However, not all inter- 



acting systems show enhanced star-forming activ- 
ity. The effect of orbital parameters and internal 
structure of the galaxies play an important part in 
the details of gas dynamics leading to star forma- 



tion (see, e.g., Cox et al. (20081; Mihos & Hern- 
|quist| (1996)). The detailed relationship between 
galaxy interactions and star formation is not well 
understood. However, in general it is found that 
the star formation rate (SFR) is increased in merg- 



ing spiral f 


galaxies (e.g., Kennicutt & Keel 


(1984); 


Keel et al. 


(1985 


); 


Bushouse (|1986 ) 


Hummel et al. 


(1990) 


, Darg et al. 


(2010b|). Therefore, we use the 



relative star formation rates in our observed sam- 
ple as an indicator of ongoing galaxy interactions. 

Finally, AGN activity may be triggered by 
galaxy interactions. This is because a huge 
amount of gas must be funneled to the nuclear 
regions to fuel the AGN's activity. There are vari- 
ous physical mechanisms by which gas loses angu- 
lar momentum and falls to the center of a galaxy, 
including gravitational torques, viscous torques, 
and hydrodynamical torques or shocks. Galaxy 
interactions and mergers provide efficient mecha- 
nisms for funneling gas to the nucleus of a galaxy 
(Mihos & Hernquist 1996). There are links be- 



tween galaxy interactions and the onset of nuclear 



activity in galaxies in both observations (Sanders 
et al.|1988 ) and simulations ( |Springel et al.|2005 1 



Many active galaxies show signs of interaction and 
merger such as disturbed morphology, particularly 
in HI (Kuo et al. 2008). Seyfert nuclei are found 



preferentially in interacting galaxies ( 


Kennicutt 


& Keel 


1984 


Laurikainen & Salo 


1995 


). The 



strongest links between interactions and AGN are 
found in the systems with the highest luminosity 
( |Bahcall et aT^1997 ). 

Using these three factors, we calculate an inter- 
action index, described below. The index includes 
only those galaxies within ± 1 Mpc of the most 
massive galaxy in each group. For galaxies with 
no distance measurement, the distance of the most 
massive galaxy in the group was assumed. The 
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index is applied to the sample of galaxy groups 
studied for this paper as well as the M81 and NGC 
2403 groups and the Local Group. Note that the 
numbers resulting from these calculations are arbi- 
trary; in order to use the index to rank our obser- 
vation sample, we have normalized each variable 
to the maximum for our sample. A robust statistic 
would need to include a normalization of the in- 
dices to some standard and appropriate weighting 
of each variable. Fortunately, for the purposes of 
this study it is only necessary to rank the groups 
relative to one another. 

We define the evidence for recent or ongoing 
interactions as follows: 



N 



Tidal 



SFR 



Nagn 



max(Nxidai) max(SFR) max(NAGN) ' 

(1) 

Here, the first factor describes the degree of dis- 
turbed morphology; the second factor the star 
formation rates, and the third factor accounts for 
AGN activity. N Tidal is the number of group 
galaxies exhibiting tidal tails or bridges. SFR is 
the average star formation rate of group galax- 
ies, calculated following Kennicutt (1998), us- 
ing Galactic extinction-corrected Ha luminosi- 
ties 11HUGS study of |Kennicutt et al.| fl2008| ). 
For galaxies with no Ha luminosity measurement 
in the 11HUGS survey ( |Kennicutt et al.||2008| ), 
the 11HUGS survey detection limit was assumed. 
A^agn is the number of group galaxies classified 
as any type of AGN in NEEQ 

We applied Equation (1) to the five galaxy 
groups. The results are tabulated in Tables [8] and 
[9j To avoid giving any one factor an artificially 
high or low weight, we assign a rank from 1-6 
for each variable, and add the ranks for the final 
ranking order for our sample. Since this ranking 
scheme is, at best, an ad-hoc indicator of ongoing 
interactions, we also tested alternate versions of 
Equation (1). We tested each permutation of two 
of the three factors, as well as each factor individ- 
ually, and each version that we tested yielded a 
slightly different result. Since none of the factors 



2 NASA/IPAC Extragalactic Database (NED) which is 
operated by the Jet Propulsion Laboratory, Cali- 
fornia Institute of Technology, under contract with 
the National Aeronautics and Space Administration, 
|http://nedwww.ipac.caltech.edu"7| 



is an iron-clad indicator of ongoing interactions in 
and of itself, we concluded that it is best to use 
all available evidence in ranking the galaxy groups 
with respect to one another. 

The M81 group has the highest interaction in- 
dex, while the NGC 45 group ranks lowest. Based 
upon these rankings, we expect that the M81 
group should contain the most HI clouds, followed 
by NGC 672, CVn I, NGC 2403, and NGC 45. 
To place these calculations in context, we have 
also calculated the interaction index for the Local 
Group. The Local Group measurements are in- 
cluded in boldface in Tables [8] and [9l We consider 
the Local Group to consist of the Milky Way, M33, 
and M31 since Local Group dwarf galaxies would 
not be visible from the distances of the groups we 
observed. We calculate the M31 and M33 star 
formation rate from the 11HUGS survey Ha lumi- 
nosity, and assume the widely-accepted star for- 
mation rate of 1 M yr~ 1 for the Milky Way. The 
Local Group ranks second to the M81 group in 
interaction evidence. 

We wish to point out a few notable areas where 
our ranking scheme is lacking. Firstly, we have 
not attempted to put a quantitative timescale on 
the index. This is because the factors we take into 
account have different (and largely unknown) du- 
rations, especially in relation to interactions, and 
combining these different timescales is a daunting 
problem. Secondly, we have no method in place for 
predicting the quantity of HI clouds that ought to 
be present or their properties based on the inter- 
action indices. Rather, we predict whether a given 
galaxy group is more likely than another to host 
HI clouds at all. Finally, we have not considered 
possible cross-correlation of the three factors. 

3. Observations 

3.1. Observations and Data Reduction 

Observations of the M81 and NGC 2403 groups 



are summarized 


in 


Chynoweth et al. ( 


2008), 


Chynoweth et al. 


( 


20091, and 


Chynoweth et al. 



([2011). We observed the CVn I, NGC 672, and 
NGC 45 groups in 32 sessions between January 
and August 2009. We combined the observing ses- 
sions into 9 maps encompassing the group galax- 
ies. Maps were made by moving the telescope 
in declination and sampling every 3' at an inte- 
gration time of 2-3 seconds per sample. Strips 
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of constant declination were spaced by 3'. Corre- 
sponding maps were made by moving the telescope 
in right ascension to form a 'basket weave' pat- 
tern over the region. Total integration time was 
approximately 58 hours for CVn I, 47 hours for 
NGC 672, and 41 hours for NGC 45. Tables [10 
through [l2] give a summary of the observations, 
including the RMS noise figures. 

The GBT data were reduced in the standard 
manner using the GBTIDL and AIPS^] data re- 
duction packages. 

In order to match our velocity resolution to 
the expected linewidths of HI clouds in the group, 
spectra were smoothed to a channel spacing of 24.4 
kHz, corresponding to a velocity resolution of 5.2 
km s _1 . A reference spectrum for each of the ob- 
servation sessions was made using an observation 
of an emission-free region, usually from the map 
edges. The reference spectrum was then used to 
perform a (signal-reference)/reference calibration 
of each pixel. The calibrated spectra were scaled 
by the system temperature, corrected for atmo- 
spheric opacity and GBT efficiency. We adopted 



the GBT efficiency equation (1) from Langston & 



Turner ( 2007 ) with a zenith atmospheric opacity 



T = 0.009. 

The frequency range observed was relatively 
free of RFI, with less than 0.5% of all spectra ad- 
versely affected. The spectra exhibiting RFI were 
identified by tabulating the RMS noise level in 
channels free of neutral hydrogen emission. Spec- 
tra that showed high RMS noise across many chan- 
nels were flagged and removed. Observations were 
gridded using the AIPS task SDIMG, which also 
averages polarizations. After amplitude calibra- 
tion and gridding, a lst-order polynomial was fit 
to line-free regions of the spectra and subtracted 
from the gridded spectra using the AIPS task 
IMLIN. Only channels in high negative velocity 
ranges, where few or no sources are expected, were 
used for the fit. This simple baseline fit was ex- 
trapolated to the positive velocity range, where 
galaxies make a baseline fit unreliable. This base- 
line velocity range yielded a flat baseline for ar- 
eas free of strong continuum radio sources, which 
is true of the majority of the region. The effec- 



3 Developed by the National Radio Astronomy Obser- 
vatory; documentation at http://gbtidl.sourceforge.net 
http:/ /www. aoc.nrao.edu/aips 



tive angular resolution, determined from maps of 
3C286, is 9.15'± 0.05'. To convert to units of flux 
density, we observed the calibration source 3C286, 
whose flux density is 14.57 ± 0.94 Jy at 1.418 GHz 



( Ott et al. 1994 ). The calibration from K to Jy was 



derived by mapping 3C286 in the same way that 
the HI maps were produced. After all corrections 
for the GBT efficiency and the mapping process, 
the scale factor from K/Beam to Jy/Beam images 
is 0.43 ± 0.03. Error estimates are difficult to ob- 
tain, given the wide range in RMS noise values 
over the maps. However, the dominant error con- 
tribution is the approximately 7% uncertainty in 
the absolute calibration for these observations. 

The mass detection thresholds were calculated 
assuming that clouds would be unresolved in the 
GBT beam, using the relation: 



= 2.36x10' 



/ D 



AV 



\Mrpc J \3y J V^ms^ 1 

(2) 

where D is the distance in Mpc, a s is the RMS 
noise in one channel, AV is the channel width, 
and N is the number of channels required for a 
secure detection. We required a cloud candidate 
to be visible in at least N — 2 channels, therefore 
the lowest velocity width that could be detected 
was 10 km s _1 . 

Mass detection thresholds are listed in Tables 
10] through [l2] HI cloud candidates were found by 



visual inspection of each channel of the baseline 
subtracted spectral line cube. Cloud candidates 
are spatially distinct from group galaxies and tidal 
streams and coherent in velocity over more than 
one channel. For each HI cloud candidate, we pro- 
duced spectra (intensity versus velocity) using the 
AIPS task ISPEC. 

4. Completeness of Study 
4.1. Velocity 

Our observations cover more than ± 700 km 
s _1 with respect to the systemic velocity of any 
group galaxy. In principle, this allows us to detect 
clouds associated with any of the mechanisms that 
create HI clouds, all of which predict cloud s to be 
within 300 km s _1 of the galaxy's velocity (Miller 
et^[2009l ). 
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4.2. Position 

Studies of the Milky Way and M31 find clouds 
anywhere within 1-50 kpc of a galaxy center 



( Wakker & van Woerden 1997 Thilker et al. 2004 1 



N-body simulations of galaxy interactions predict 
that clouds and large-scale structure will be found 
as close as 3 kpc and as far as 750 kpc from galax- 
ies (Miller et al. 2009). In order to balance the 



requirement of a good mass sensitivity and a large 
angular area in our observations, we designed our 
observations to cover a projected radius of at least 
40-50 kpc around major group galaxies. The map 
sizes are adequate to detect both interaction rem- 
nants and DM halo-associated clouds. The an- 
gular resolution of our observations ranges from 
10-20 kpc. Since galactic fountain HI clouds are 
found at a distance of < 10 kpc from the galac- 
tic disk, we may not be able to resolve galactic 
fountain clouds, depending on their velocity and 
orientation. 

4.3. Mass 

The 5ct detection threshold for HI masses 
ranged from 9.7 x 10 5 to 12.7 x 10 6 M , assuming 
a linewidth of 10 km s -1 and using the distance of 
the most massive group galaxy. With these mass 
detection thresholds, we would be able to detect 
analogs of the most massive HI clouds that have 
been detected surrounding nearby galaxies, as well 
as large Milky Way objects such as Complexes C 
and H and the Magellanic Stream. These clouds 
are the most likely to originate in galaxy inter- 
actions, cold accretion, and dark matter halos, 



rather than from supernovae (Miller et al. 20091, 



so our study is biased towards HI clouds with an 
extragalactic origin. To estimate the complete- 
ness of the survey in terms of mass, we compiled 
the published values for all known extraplanar HI 
clouds surrounding nearby galaxies and the Milky 
Way. If this mass spectrum is representative of 
the general extragalactic HI cloud mass function, 
we can detect anywhere from 12% to 35% of the 
analogous clouds surrounding the galaxies in our 
sample. Table [l3| lists the average mass detection 
threshold and mass completeness for each galaxy 
group. Figure [2] shows the mass function of known 
HI clouds compared to our observational limits for 
each galaxy group. Note that the mass function 
for HI clouds is not well-constrained, particularly 



for the Milky Way, since the properties of Milky 
Way HI clouds beyond the stellar halo are not 
known. We should be able to detect cold accre- 
tion clouds, which are predicted to have HI m asses 
of 10 6 to 10 7 M flKeres fc Hernquist||2009 ). 



5. Predictions of Numerical Simulations 

In order to determine properties expected for 
gas-embedded dark-matter halos, we analyzed sev- 
eral N-body simulations. We conducted two dark 
matter-only cosmological N-body simulations for 
comparison. One simulation uses initial conditions 
from WMAP3, and one from WMAP5. We simu- 
lated a 50 3 (Mpc/h) 3 volume of the universe from 
z=149 to z=0 using 256 3 particles. At redshift 
zero, we selected a volume 10 Mpc/h on a side 
to resimulate with 512 3 particles with a 'zoom' 
technique aimed to preserve the tidal field out- 
side the higher resolution volume. This smaller 
volume was selected to host a ~ 10 12 M Q halo at 
z=0, consistent with the M81 galaxy halo. We 
used a friends-of-friends algorithm with a linking 
length of b = 0.2 to identify dark matter halos 
(IDavis et al.||1985b, and used a SUBFIND tech- 



nique to determine the bound subhalos (Springel 
etaT1[200H ). 



We also use the dark matter-only "set C" simu- 



lations from Polisensky & Ricotti (2011 1 in which 
a 65 3 (Mpc/h) 3 box was simulated with 408 3 parti- 
cles and WMAP3 cosmological parameters. Nine 
dark matter halos with virial masses at z=0 of 
1.6-2.5 x 10 12 M were resimulated using a zoom 
technique with mass resolution similar to the sim- 
ulations described above. The AHF halo finder 



( Knollmann & Knebe 2009 ) was used to determine 



the bound subhalos. 

In order to connect the numerical simulations 
with the properties of HI observations, we require 
an estimate of the HI to dark matter mass frac- 
tion. The distribution of HI gas compared to the 
dark matter halo distribution is uncertain. For the 
purpose of comparing our HI observations to the 
simulated HI clouds, we assume that HI gas fol- 
lows the dark matter distribution. We include in 



our comparison the prescription of Gnedin ( 2000 1 , 



which includes a mass cutoff such that a dark mat- 
ter halo with mass less than 2 x 10 s M has no 
associated HI gas. We take as the scaling factor 
from dark matter to HI gas mass the|Zwaan et al.| 
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(2003) results, summarized in Fukugita & Peebles 



(2004). Their values yield a HI to dark matter 



mass fraction of 0.0018 ± 0.0003. 

After scaling the simulated dark matter to HI, 
we apply the 5er HI mass threshold, velocity lim- 
its, and angular size limits of our survey. For each 
simulation, we calculate the probability that de- 
tectable sub-halos reside around the major dark- 
matter halos. This probability depends on the 
mass sensitivity of each observation, ranging from 
38% of major dark-matter halos hosting detectable 
sub-halos at the mass sensitivity of the NGC 45 
group observation, to 97% at the M81 and NGC 
2403 groups mass sensitivity. For each major 
dark matter halo with sub-halos that would be 
detectable as HI clouds in our observations, we 
find the properties of those sub-halos. We repeat 
this process rotating each simulation box three 
ways. In the simulations, we find between and 
16 clouds with an average of 3 clouds of sufficient 
mass to be detected in our observations. Figure 
[3] shows the distribution of number of detectable 
simulated HI clouds per major dark-matter halo. 
It is important to note that in the NGC 672 
and NGC 45 groups, the probability of the ma- 
jor halo hosting a dark-matter associated HI cloud 
detectable by our observations is < 50%. Figure 
[4] shows the position- velocity distributions of the 
simulated HI clouds for each galaxy group, and 
Figure [5] shows the simulated HI mass functions. 
The simulated clouds have a wide distribution of 
position and velocity separations from the central 
galaxy, ranging out to the very edge of our obser- 
vations. 

Table[T4|presents the detailed predictions of the 
simulations. It is worth noting that the predicted 
properties (positions and velocities) of the clouds 
are remarkably similar, regardless of the varied 
mass sensitivities of the observations. The sensi- 
tivity of the observation appears to affect only the 
number of clouds we should detect, and not their 
properties. That is, the positions and velocities 
of the simulated HI clouds do not correlate with 
their masses. Therefore, we would not expect that 
more sensitive observations would find HI clouds 
with similar position and velocity distributions to 
those observed in the M81 group HI clouds. No HI 
clouds were detected in the corresponding region 
of our observations of the CVn I, NGC 672, and 
NGC 45 groups. 



The above analysis does not account for numer- 
ical destruction effects in the simulations. Due 
to the gravitational softening lengths employed in 
N-body simulations subhalos orbiting within the 
halo of a larger galaxy experience extraneous tidal 
fields than can artificially disrupt subhalos. These 
effects are dominant for poorly resolved subhalos 
in the inner halo region that have short dynami- 
cal times. A method for estimating the amount of 
numerical destruction in simulations is presented 



in (Polisensky & Ricotti 20111 and applied to 



their 9 set C halos the method estimates that be- 
tween 2 to 6 subhalos with dark matter masses 
> 2 x 1O 8 M0 were numerically destroyed within 
50 kpc of each halo center. This would increase the 
predicted number of detectable subhalos and make 
the disparity with the observations more striking. 

However, our dark matter-only simulations do 
not include the effects of baryons on the dark mat- 
ter subhalo abundances. Dark matter subhalos 
whose orbits intersect a baryonic disk can be de- 



stroyed by disk shocking (Ostriker et al. 19721 



Calculations show a reduction in the dark matter 
subhalo abundances by factors of 2 — 3 for masses 
10 9 — 10 7 M Q in the inner ~ 30 kpc (~ 5 disk scale- 
lengths) of a Milky Way-sized galaxy with a disk 
( D'Onghia et al.||2010 ). Disk shocking is likely to 
offset the gains in subhalo numbers from consider- 
ations of numerical destruction for group galaxies 
with disk structures. 



In Chynoweth et al. 



(20111, we presented ob- 
servations of the M81 and NGC 2403 groups. We 
found 9 HI clouds associated with the M81 group 
so in this case a statistical comparison was possi- 
ble. Chynoweth et al. (2011) presents a detailed 
description of the clouds. Table [15] summarizes 
the properties of the clouds in the M81 group for 
reference. To summarize, the detected HI clouds 
have masses of 10 6 ~ 8 M Q , and have a mean veloc- 
ity deviation of 125 km s _1 and mean projected 
distance of 110 kpc from M81. In contrast, analy- 
sis of numerical simulations (in the same manner 
described in this section) predicts that we should 
detect 40-50 HI clouds embedded within dark mat- 
ter minihalos. These simulated clouds would have 
masses of 10 6-8 M©, average velocity deviation 
from the simulated "M81" galaxy of 250 km s _1 , 
and mean projected separation of 170 kpc. Al- 
though the detected HI cloud mean masses match 
the simulated mean masses, the observed and sim- 
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ulated mass functions are discrepant. A KS-test of 
the two mass functions yields only a 0.3% proba- 
bility that the mass functions are statistically sim- 
ilar. Additionally, the number count and distri- 
butions of positions and velocities of detected HI 
clouds are statistically inconsistent with the distri- 
butions expected for dark matter halos with em- 



bedded HI based on the simulations (Chynoweth 
et al.poTT ). For the groups observed for this pa- 
per, no HI clouds were detected so the discrepancy 
between predictions and observations is further in- 
creased. 

6. Results, Conclusions and Future Work 

We have observed 5 nearby galaxy groups with 
the GBT in an effort to search for HI clouds as- 
sociated with either galaxy interactions or dark- 
matter halos. In Chynoweth et al. ( |2011[ ), we pre- 
sented observations of the M81 and NGC 2403 
groups. We found 9 HI clouds associated with 
the M81 group, and no clouds in the NGC 2403 
group. In this study, we observed the CVn I, NGC 
672 and NGC 45 groups. The 5 observed groups 
have varying levels of galaxy interactions, which 
we quantify using an interaction index. 

Figures [6] and J7] show the CVn I Group obser- 
vations. Figures |8pmd [9] show the NGC 672 group 
observations. Finally, Figure [10] shows the NGC 
45 group observational results. In these HI column 
density (moment 0) maps, we show all HI over the 
velocity ranges spanning each galaxy group, with 
galaxies labeled. Note that in some maps, there 
appears to be emission not associated with known 
galaxies. All apparent extra-galactic emission was 
found not to have coherent velocity structure, so 
this emission is likely due to noise such as unde- 
tected RFI or edge effects from the gridding func- 
tion. 

We detected no HI clouds in the CVn I, NGC 
672 or NGC 45 groups above our mass detection 
thresholds and within ± 700 km s _1 and 50 kpc 
of group galaxies. It is important to note that the 
NGC 672 and NGC 45 groups are the most distant 
of our observed sample, and therefore our observa- 
tions are less sensitive for these groups (see Figure 
|2]). More sensitive observations of all these groups 
would be necessary in order to better constrain 
their HI properties. 

The observed population of HI clouds and 



HVCs around the Milky Way and other nearby 
galaxies, including those we detect in the M81 
group, arc likely created through a combination of 
galactic fountains, galaxy interactions, dark mat- 
ter substructure, and cold accretion. The focus of 
this paper is to address the role of galaxy inter- 
actions and gas-embedded dark matter halos as 
originating mechanisms for clouds with HI masses 
greater than ~ 10 6 M . We find the role of each 
mechanism to be as follows. 

The observed M HI > 10 6 M Q HI cloud popu- 
lation does not match the predictions of ACDM 
cosmological models. We have analyzed a suite of 
cosmological N-body simulations in order to pre- 
dict the number and properties of such clouds. 
Our simulations predict that the number of de- 
tected HI clouds in each group should have ranged 
between 0-16, with an average of 3 HI clouds de- 
tected per group with our observations. The HI 
clouds that were detected in the M81 group are 
not statistically coincident in position and veloc- 



ity with the simulated clouds (Chynoweth et al. 



2011). Unfortunately, given the mass sensitivity 



of our observations in the NGC 672 and NGC 45 
galaxy groups, we cannot say with confidence that 
dark-matter associated HI clouds ought to have 
been detected. It is possible that the rest of the 
groups in our sample are all among those major 
halos in the simulations that hosted no detectable 
HI clouds-our results would become more robust 
with more observed galaxy groups and/or more 
sensitive observations of this sample. 

We find HI clouds only in the M81 galaxy 
group, which is currently undergoing strong in- 
teractions and ranks highest in our "interaction 
evidence" ranking system. Our results therefore 
indicate that the majority of Mhi > 10 6 M ex- 
tragalactic HI clouds are generated through tidal 
stripping caused by galaxy interactions. By anal- 
ogy, the Milky Way HVC population is likely gen- 
erated through a combination of tidal stripping 
and galactic fountains (galactic fountains in the 
galaxy groups being un-detectable in our observa- 
tions). 

Ideally, we would also address cold accretion as 
a possible source of HI clouds. Currently, there 
is no published catalog of the properties of cold 
accretion-generated HI clouds, and therefore we 
cannot make a definitive statement about the role 
of cold accretion. If HI clouds tracing cold accre- 
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tion exist in our observed position, velocity, and 
mass regimes, as the simulation of Kere s fe Hern-| 
quist ( 2009 ) predict they ought (see Table[l]), they 
should have been observed in all groups instead of 
only one group. If the HI clouds we have detected 
are tracers of cold accretion, this result indicates 
that the presence or lack of cold accretion-related 
HI clouds is strongly environment-dependent, re- 
quiring a dense environment and/or the presence 
of current, strong interactions such as those in the 
M81 group. Tentatively, we state that the ob- 
served HI clouds are not likely to be a product 
of cold accretion. Certainly, more modeling and 
simulations of cold accretion, especially in the very 
cold HI regime, are needed to properly address this 
question. 

Future work towards determining the origin of 
HI clouds in galaxy groups should include high- 
resolution HI observations with the EVLA to de- 
termine their substructure; UV and optical obser- 
vations to search for a stellar component; and com- 
parison with simulations of the M81 Group that 
include all the relevant gas physics and allow for 
cold accretion. These studies would place stronger 
constraints on the origins of HI clouds in galaxy 
groups, and determine whether such clouds are 
analogs to the Milky Way HVCs. Finally, these 
future observations will place stronger constraints 
on cosmological models of galaxy formation. 
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Fig. 1. — Map of observed galaxy group locations. The gray wedge indicates the declination range not 
observable from Green Bank. 
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Fig. 2. — Mass distribution of known HI clouds associated with the Milky Way and nearby galaxies M31, 
M33, M1Q1, M83, and NGC 2997. Our mass detection thresholds are indicated with vertical lines. Clouds 
to the right of the lines would be detectable in our observations. If this population of clouds is representative 
of galaxy groups, our detection limit allows us to detect the most massive 11-34% of HI clouds in the groups. 
Data: 



et al. (2006); Simon et al. 



Thilker et al.| {2004} ; |Grossi et aL| p008| ; |van der Hulst fc Sancisi| (|1988|; |Wa kker et al.| p007 ) 
(|2006[); |van Woerden et al.| (|1999[); |Putm an et al.| ( |2003| ) 



Braun & Thilker (2004); Miller et al. (2009); Hess et al. (2009). 



Lockman et al. (20081; 



Thorn 



Table 1: Expected Properties of HI Clouds 



HI Cloud Origin 


Distance 
(kpc) 


Deviation Velocity 
(km s" 1 ) 


Size 
(kpc) 


Mass 
M 


Galactic Fountain a 
Galaxy Interactions 3, 
Dark Matter Halos a 
Cold Accretion* 3 


< 10 
5-100 
150-750 
30-60 


< 150 

< 300 

< 300 


< 1 
~ 1 
1-10 


< 10 4 

10 5 - 10 6 

10 6 - 10 8 
10 6 -10 7 







"From 
'From 



Miller et af7| ||2009|) 



Keres & Hernquist I 2009 1 



16 



Number of Simulated HI Clouds per Holo 
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Number of Simulated HI Clouds per Halo 




Nciouds Per halo 
Number of Simulated HI Clouds per Halo 

a5 r— , , , . , , , , , , , , , , . r- 



| 

'in 



0.4 



0.3 



0.2 



0.1 



0.0 



CVn 



2 4 6 8 

N Clouds P er hal ° 

Number of Simulated HI Clouds per Halo 
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0.6 

0.5 

0.4 

0.3 

0.2 

0.1 
0.0 



NGC 672 



2 4 6 

Nciouds Per halo 



Fig. 3. — Number distribution function of simulated clouds over all simulations described in Section [5] 
Diagrams include all dark matter halos above the mass detection threshold for the galaxy group (assuming 
a HI to dark matter mass fraction of 0.0018 ± 0.0003 as described in Section [5]), that are within ± 700 km 
s _1 and ± 50 kpc of the most massive galaxy in each galaxy group. The distributions plotted here show the 
number of subhalos per major halo. The distribution functions are normalized so that their sums are unity. 



17 




io - 



E^i . u^,^ i i 1 I I 0.0 

D 100 200 300 400 500 600 700 

IV-V(Center)! (km s'') 

Fig. 4. — Position-velocity distributions of simulated clouds. Distributions include all simulated dark matter 
halos above the mass detection threshold for the galaxy group (assuming a HI to dark matter mass fraction 
of 0.0018 ± 0.0003 as described in Section [HJ, that are within ± 700 km s _1 and ± 50 kpc of the most 
massive galaxy in each galaxy group. The greyscale for each plot runs from zero to the maximum number 
of subhalos detected in that simulated galaxy group, normalized so that the maximum is unity. 
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Fig. 5. — HI mass functions of simulated clouds. Distributions include all simulated dark matter halos above 
the mass detection threshold for the galaxy group (assuming a HI to dark matter mass fraction of 0.0018 ± 
0.0003 as described in Section pL that are within ± 700 km s _1 and ± 50 kpc of the most massive galaxy 
in each galaxy group. The mass distribution functions are normalized so that their sums are unity. 
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Fig. 6. — CVn I Group HI column density I: Group galaxies and background galaxies are labeled. 
All flux above la is included. Grayscale: 0-150 Jy beam -1 x km s _1 . Contours: 2 Jy beam -1 x km s _1 
x (-3,3,5,10,25,50,100,250) 
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Fig. 7. — CVn I group column density maps II and III. Left: HI column density map of M94. 



Grayscale: 0-50 Jy beam x km s 
HI column density map of NGC 4395. 
kms" 1 x (-3,3,5,10,25,50,100,250) 



Contours: 2 Jy beam -1 x kms" 1 x (-3,3,5,10,25,50,100,250. Right: 
Grayscale: 0-75 Jy beam -1 x km s _1 . Contours: 2.75 Jy beam -1 x 



Table 2: Galaxy Groups Within 7.5 Mpc 



Name" 


a 


S 


D 






(J2000) 


(J2000) 


(Mpc) 


(km s -1 ) 


M81 


10:10 


+69.1 


3.7 


127 


Maffei 


03:49 


+68.0 


3.1 


171 


Sculptor 


00:36 


-31.0 


2.8 


177 


NGC 2403 


07:43 


+66.7 


3.7 


231 


CVn I 


12:18 


+35.8 


4.0 


294 


NGC 5128 


13:25 


-41.6 


4.3 


390 


NGC 672 


01:51 


+27.8 


6.1 


458 


NGC 45 


00:13 


-23.4 


6.1 


503 



a Galaxy groups observed for this study are listed in bold. 
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Fig. 8. — NGC 672 Group HI column density I: Any emission not associated with labeled galaxies is 
likely due to noise, since it does not have coherent velocity structure. Grayscale: 0-115 Jy beam -1 x km 
s -1 . Contours: 3.25 Jy beam -1 x km s -1 x (-3,3,5,10,25,50,100,250) 
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Fig. 9.— NGC 672 Group HI column density II and III. Left: HI column density map of NGC 784. 
Grayscale: 0-25 Jy beam" 1 x kms -1 . Contours: 1.0 Jy beam" 1 x km s _1 x (-3,3,5,10,25,50,100,250. Right: 
HI column density map of UGC 1281. Grayscale: 0-25 Jy beam -1 x km s _1 . Contours: 1.0 Jy beam -1 x 
km s" 1 x (-3,3,5,10,25,50,100,250). Any emission not associated with labeled galaxies is likely due to noise, 
since it does not have coherent velocity structure. 
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Fig. 10. — NGC 45 group HI column density. Any extraplanar emission is likely due to noise, since 
it does not have coherent velocity structure. Top Left: HI column density map of NGC 45. Grayscale: 
0-100 Jy beam" 1 x km s" 1 . Contours: 0.30 Jy beam" 1 x km s" 1 x (-3,3,5,10,25,50,100,250. Top Right: 
HI column density map of NGC 24. Grayscale: 0-50 Jy beam -1 x km s _1 . Contours: 0.75 Jy beam -1 x 
km s" 1 x (-3,3,5,10,25,50,100,250). Bottom Left: HI column density map of NGC 59. Grayscale: 0-2.5 Jy 
beam" 1 x kms" 1 . Contours: 0.15 Jy beauT 1 x km s -1 x (-3,3,5,10,25,50,100,250) 
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-12.60 ± 0.03 
-12.61 ± 0.05 
-13.16 ± 0.11 
-13.36 ± 0.05 
-13.70 ± 0.16 
-12.70 ± 0.07 


Type 


SAB(s)dm 
Im 
Im 

SAB(s)dm IV 
Im: 
ImlV-V 
IAm 


> a 


lO ^ N N n n o 
CO N O Ol M lO lO 

PI M M H M H H 


o 
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3.33 
4.22 
3.13 


(J2000) 
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24:47:00 
40:44:00 
38:01:00 
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23:03:00 
44:31:00 


(J2000) 


13:34:42 
13:39:18 
13:39:54 
13:50:54 
13:54:48 
14:15:54 
14:24:42 


Name 


NGC 5238 
UGC 8638 
UGC 8651 
UGC 8760 
UGC 8833 
UGC 9128 
UGC 9240 
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7 

'a 

5^ SB, 
ao 

ao 


-11.87 ± 0.04 
-11.22 ± 0.04 
-12.36 ± 0.06 


Type 


SA(s)c 
SA(s)dm 
SA(rs)0-: 


> a 


-tf b- CN 

in o o 
m ^ co 




7.88 
6.74 
5.37 


(J2000) 


-24:57:47 
-23:10:55 
-21:26:40 


a 

(J2000) 


00:09:56.5 
00:14:04.0 
00:15:25.1 


Name 


NGC 24 
NGC 45 
NGC 59 



CO 



30 



66 



mean SFR/max(mean SFR) 


1 

0.22 
0.05 
0.06 
0.15 
0.14 


fa i 

CO H 

03 
CD fe^ 


0.47 

0.103 
0.023 
0.027 
0.071 
0.067 




_ LO 

O ^ i-< O O 




O i — 1 CN CN O O 


h 

¥ 

a 

"a 

-a 

h 
55 


lO LO 

. • o o o 
o o 


"3 

-a 

h 
55 


cm o o o 


Galaxy Group 


Local Group 

M81 

NGC 672 
CVn I 
NGC 2403 
NGC 45 
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Tabic 9: Interaction Index: Ranked 



Galaxy Group 






Mean SFR 


I 


I Ranked 


M81 


6 


5 


5 


16 


1 


Local Group 


4 


1 


7 


12 


2 


NGC 672 


4 


6 


1 


11 


3 


CVnl 


1 


7 


2 


10 


4 


NGC 2403 


1 


1 


4 


6 


5 


NGC 45 


1 


1 


3 


5 


6 



Table 10: CVn I Group Observations Summary 
Area: 



NGC 4244 a range (J2000): 


12:05:00.0 - 


12:25:00.0 


NGC 4244 S range (J2000): 


35:00:00.0 - 


39:00:00.0 


M94 a range (J2000): 


12:46:00.0 - 


12:56:00.0 


M94 S range (J2000): 


40:00:00.0 - 


42:00:00.0 


NGC 4395 a range (J2000): 


12:21:00.0 - 


12:31:00.0 


NGC 4395 S range (J2000): 


32:30:00.0 - 


34:30:00.0 


Observations: 






Center Frequency (MHz): 


1420.4 




Bandwidth (MHz): 


12.5 




Velocity Range (Heliocentric, km s -1 )): 


-1315 - 1320 


Channel Width (kHz): 


24.4 




Velocity Resolution (km s -1 -* 


5.2 




Integration time (hours): 






NGC 4244: 


42.2 




M94: 


9.8 




NGC 4395: 


5.7 




Typical RMS noise (mK/channel): 






NGC 4244: 


26.3 




M94: 


21.2 




NGC 4395: 


32.6 




Sensitivity to HI (ler): 






M HI (x 10 5 M ) 






NGC 4244: 


2.40 




M94: 


4.16 




NGC 4395: 


4.27 




Beam size (kpc): 






NGC 4244: 


12.7 




M94: 


13.9 




NGC 4395: 


11.3 
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Table 11: NGC 672 Group Observations Summary 
Area: 

NGC 672 a range (J2000): 01:38:00.0 - 01:52:00.0 

NGC 672 S range (J2000): 28:00:00.0 - 26:00:00.0 

NGC 784 a range (J2000): 01:58:00.0 - 02:04:00.0 

NGC 784 S range (J2000): 29:30:00.0 - 28:10:00.0 

UGC 1281 a range (J2000): 01:47:00.0 - 01:15:00.0 

UGC 1281 S range (J2000): 33:10:00.0 - 32:00:00.0 

Observations: 

Center Frequency (MHz): 1418.4 

Bandwidth (MHz): 12.5 

Velocity Range (Heliocentric, km s -1 )): -890 - 1740 

Channel Width (kHz): 24.4 

Velocity Resolution (km s -1 ' 5.2 

Integration time (hours): 

NGC 672: 36.4 

NGC 784: 5.3 

UGC 1281: 5.7 

Typical RMS noise (mK/channel): 

NGC 672: 5.92 

NGC 784: 12.3 

UGC 1281: 12.8 

Sensitivity to HI (ler): 

M H i (x 10 5 Mq) 

NGC 672: 1.99 

NGC 784: 1.88 

UGC 1281: 2.41 

Beam size (kpc): 

NGC 672: 19.8 

NGC 784: 12.2 

UGC 1281: 13.6 
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Table 12: NGC 45 Group Observations Summary 



A TPfi • 
ill cel. 






atpp ran OA ( T9000V 


00:11:00.0 - 


00:17:00.0 


inLtO 4o o range (^jzuuuj 


-22:30:00.0 


- -23:50:00.0 


lUPP OA n ranra ( 79000 V 


00:07:00.0 - 


00:13:00.0 


WPP 9/1 A rano-A f 79000V 


-24:20:00.0 


- -25:40:00.0 


NOP PiQ rv ranp-e f 79000V 


00:13:00.0 - 


00:18:00.0 


AJPP X rono-A ( 79000V 


-20:50:00.0 


- -22:10:00.0 


Observations: 








1420.4 




RanHwiHtri (MH7V 


12.5 




a 1 \^\t Rfincrp f77plirir'0'n1~T - ir' km Q ~^V 
v (-■lULji L^y j. vciiigjV^ iaa^iiul-i^iilji . iyiii 11. 


-1315 - 1320 


PViannol WiH+Vi /"H7r?V 


24.4 




Velocity Resolution (km s 1 -' 


5.2 




Integration time (hours): 








11.3 




t\tpp OA- 


12.9 




inLtO oy 


16.5 




_L V JJl^Oil J. V,lV_lk_J IIUIdL I 111 J. V / ^IIOjIIIUjI I . 






NGC 45: 


13.9 




NGC 24: 


10.2 




NGC 59: 


10.9 




Sensitivity to HI (la): 






M HI (x 10 5 M Q ) 






NGC 45: 


4.72 




NGC 24: 


4.72 




NGC 59: 


2.19 




Beam size (kpc): 






NGC 45: 


21.3 




NGC 24: 


18.2 




NGC 59: 


14.5 





Table 13: Mass Detection Thresholds 



Galaxy Group 


Average Detection Threshold 
(5a, x 10 6 M ) 


Fraction of Detectable HI Clouds 


M81 


0.97 


0.34 


NGC 2403 


1.03 


0.32 


CVn I 


5.23 


0.23 


NGC 672 


6.43 


0.21 


NGC 45 


12.7 


0.11 
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Table 14: Numerical Simulation Predictions' 



Galaxy Group 


Probability of 


< N c i ou a s > 


< M H i > 


^ Vcloud Vgalaxy ^* 


Rcloud ^tgalaxy 




detectable HI clouds 




(x 10 6 M ) 


(km s- 1 ) 


(kpc) 


M81 


0.97 


6 


7.4 


172 


37 


NGC 2403 


0.97 


5 


7.0 


209 


36 


CVn I 


0.67 


2 


19.8 


201 


35 


NGC 672 


0.50 


2 


11.4 


251 


35 


NGC 45 


0.38 


1 


39.3 


242 
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a All values calculated within ± 50 kpc and ± 700 km s 1 of central halo, and averaged over all numerical simulations. Properties 
of HI clouds are calculated for the major dark-matter halos that host detectable HI clouds. 



Table 15: M81 Filament HI Cloud Properties 



Cloud a 


Coordinate Designation 


ADmsi 
(kpc) 


Tpeak 

(K) 


AVmsi 
(km s- 1 ) 


a v 
(km s _1 ) 


M HI ( 3 . 63Mpc ) 

x 1O 6 M 


C08-1 


GBC J095007.2+695556 


65 


0.11 


202 


50 


14.7 


C08-2 


GBC J100250. 7+681949 


63 


0.10 


68 


55 


22.5 


C08-3 


GBC J095244.9+681250 


57 


0.12 


48 


82 


26.7 


C08-4 


GBC J100145.0+691631 


38 


0.30 


108 


28 


83.7 


C08-5 


GBC J095506.3+692205 


21 


0.07 


317 


36 


6.9 


C10-1 


GBC J092635.8+702850 


186 


0.13 


144 


57 


3.6 


C10-2 


GBC J101926.3+675222 


161 


0.20 


72 


26 


3.4 


C10-3 


GBC J091952. 7+680937 


218 


0.17 


74 


42 


5.1 


C10-4 


GBC J102239. 1+684057 


159 


0.30 


103 


68 


12.0 


Mean Values 
Simulated Mean Values b 




108 
170 


0.17 


126 
250 


49 


19.8 
10.8 



a C08 indicates clouds from Chy nowe th et ah] |2008 ) , C09 indicates clouds from [Chynoweth et al.| | |2009fl , and C10 indicates 

cloud s from | Chynoweth et al. ( |20 lip . 
''From [Chynoweth et aL (2011 1: calculated in a region matching the large size (± 250 kpc) and velocity range (± 900 km s — 1 ) 

of the M81 Filament observation. 
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